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Abstract 

We consider neutrino oscillations and neutrinoless double beta decay in a 
five dimensional standard model with warped geometry. Although the see- 
saw mechanism in its simplest form cannot be implemented because of the 
warped geometry, the bulk standard model neutrinos can acquire the desired 
(Majorana) masses from dimension five interactions. We discuss how large 
mixings can arise, why the large mixing angle MSW solution for solar neutri- 
nos is favored, and provide estimates for the mixing angle U e 3- Implications 
for neutrinoless double beta decay are also discussed. 
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1 Introduction 



Extra dimensions are an interesting possibility for physics beyond the standard 
model (SM). They offer an explanation for the weakness of gravity by the large vol- 
ume of the compactified space or by localization of the graviton at a space-time 
boundary Gravity becomes strong at the TeV scale, and the large hierarchy 

between the scale of gravity and the weak scale is eliminated. 

We take the fifth dimension to be an S1/Z2 orbifold with a negative bulk cosmo- 
logical constant, bordered by two 3-branes with opposite tensions and separated by 
distance R. Einstein's equations are satisfied by the non-factorizable metric 

ds 2 = e- 2a{y) r]^dx^dx u + dy 2 , a(y) = k\y\ (1.1) 

which describes a slice of AdSs. The 4-dimensional metric is rj^ = diag(— 1, 1, 1, 1), 
k is the AdS curvature related to the bulk cosmological constant and brane tensions, 
and y denotes the fifth coordinate. The AdS curvature and the 5d Planck mass M 5 
are both assumed to be of order Mpi ~ 10 19 GeV. The AdS warp factor Q = e~ nky 
generates an exponential hierarchy of energy scales. If the brane separation is kR ~ 
11, the scale at the negative tension brane, located at y = ttR, is of TeV-size, while 
the scale at the brane at y = is of order Mpj. At the TeV-brane gravity is weak 
because the zero mode corresponding to the 4D graviton is localized at the positive 
tension brane (Planck-brane). 

Models with localized gravity open up attractive possibilities for flavor physics. If 
the SM fermions reside in the 5-dimensional bulk, the hierarchy of quark and lepton 
masses can be interpreted in a geometrical way @,||. Different flavors are localized 
at different positions in the extra dimension or, more precisely, have different wave 
functions. The fermion masses are in direct proportion to the overlap of their wave 
functions with the Higgs field ||. Moreover, bulk fermions reduce the impact of 
non-renormalizable operators which, for instance, induce rapid proton decay and 
large neutrino masses @-|U- 

In the proposed framework it is natural to take the SM gauge bosons as bulk 
fields as well to maintain 5D gauge invariance. The Higgs field has to be localized at 
the TeV-brane. Otherwise the gauge hierarchy problem reappears J7L ^| . Comparison 
with electroweak data, in particular with the weak mixing angle and gauge boson 
masses, requires the Kaluza-Klein (KK) excitations of SM particles to be heavier 
than about 10 TeV ||||]. If the fermions were confined to the TeV-brane, the KK 
scale would be even more constrained [p|[iO|.Pl 

In refs. [0,0 it was demonstrated that small Dirac neutrino masses can be 



generated by adding sterile neutrinos in the bulk. Both the solar and the atmospheric 
neutrino anomalies can be accommodated. 



3 Recently it has been advocated that the third generation of quarks should be confined to the 
TeV-brane to prevent large contributions to the electroweak p parameter [ [TT| . However, no such 
conclusion can be drawn for the third generation of leptons. 
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In this letter we investigate the alternative possibility of Majorana neutrinos. 
Because of the warped geometry, the usual see-saw mechanism cannot be imple- 
mented in a straightforward way. Namely, under the assumption that the large 5D 
Majorana mass is y independent, the lowest state of the singlet neutrino, it turns 
out, is localized close to the TeV-brane where it acquires a typical KK scale mass 
(~ a few TeV) from the warped geometry. Fortunately, in contrast to the SM, the 
warped 5D SM admits larger dimension five Majorana masses and, with a judicious 
choice of parameters these can be employed for resolving the solar and atmospheric 
neutrino anomalies. Large mixings arise if the SM neutrinos are located at similar 
positions in the extra dimension. We demonstrate that the large mixing angle MSW 
is the favored solution to the solar neutrino anomaly, while U e % can be kept below 
the observational bound. We also discuss possible experimental signatures to distin- 
guish between the scenarios of Dirac and Majorana neutrinos, such as neutrinoless 
double beta decay and /i — > ej. We also comment on the issue of proton decay. 

2 Fermions in a Warped Background 

To fix the notation let us briefly summarize the properties of fermions in a warped 
geometrical background. Since the 5D theory is non-chiral, every Weyl fermion of 
the SM has to be associated with a Dirac fermion in the bulk. The number of 
fermionic degrees of freedom is doubled. Chirality in the 4D low energy effective 
theory is restored by the orbifold boundary conditions. 

The Dirac equation for a fermion in curved space-time reads 



where E^f is the fiinfbein, 7° = (7 / \7 5 ) are the Dirac matrices in flat space, 



is the spin connection, and a' = da/dy. The index M refers to objects in 5D 
curved space, the index a to those in tangent space. Fermions have two possible 
transformation properties under the Z2 orbifold symmetry, *&(— y)± = ±'j^(y)±. 
Thus, \l/±\l/± is odd under Z2, and the Dirac mass parameter, which is also odd, can 
be parametrized as m^, = co' . The Dirac mass should therefore originate from the 
coupling to a Z 2 odd scalar field which acquires a vev. On the other hand, is 
even. Using the metric ( |1 . 1|) one obtains for the left- and right-handed components 
of the Dirac spinor [|, |T^ ] 



(2.2) 




(2.3) 



[e 2 °d^ + d 2 5 - o'd 5 - M 2 }e- 2 ^ L , R = 0, 



(2.4) 



where M 2 



c(c ± l)k 2 =f co" and ^ L , R = ±^l,r- 
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Decomposing the 5d fields as 



1 CO 

= -=X)* (n) M/n(y), (2-5) 



eq. (|2.4|) admits a zero mode solution [|, [H 



P (2-c)<7 



and a tower of KK excited states 



e 5a/2 



fn(v) = ^ [U^e°) + 6 Q (m„)F Q (^)] . (2.7) 

The order of the Bessel functions is a = |c± 1/2 1 for ^l,^- The spectrum of KK 
masses m n and the coefficients b a are determined by the boundary conditions of the 
wave functions at the branes. The normalization constants follow from 

ttR 

-3cr 



2nR 



dy e a f m (y)fn(y) = 5 mn . (2. 

-ttR 



Because of the orbifold boundary conditions, the zero mode of (^ '_) is a 
left-handed (right-handed) Weyl spinor. For c > 1/2 (c < 1/2) the zero mode of the 
fermion is localized near the boundary at y = (y = ttR), i.e. at the Planck- (TeV-) 
brane. The KK states are always localized at the TeV-brane. 



3 The See-Saw Mechanism and Warped Geome- 
try 

In the see-saw mechanism lepton number is broken by the Majorana mass of heavy 
right-handed neutrinos, Cm = M^NN. Together with a Dirac mass term Co = 
Xn^lNH, where H denotes the SM Higgs field, a Majorana mass for the SM neu- 
trinos v L is generated, 

M v = (3.9) 

M N 

Taking M v ~ 50meV (of the order of the atmospheric Am 2 ), one finds ~ 
A^ • 6 x 10 14 GeV. For 0.01 < A^ < 1 this points to an intermediate scale for the 
right-handed Majorana mass. 

In theories with TeV-scale gravity Mn is naturally bounded to be below a few 
TeV. Therefore, very small Yukawa couplings A at < 10~ 6 are required to generate 
sub-eV neutrino masses. Keeping in mind the large hierarchy of Yukawa couplings 
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for quarks and charged leptons, this result may not seem "completely unreasonable" . 
Nonetheless, a deeper understanding of the smallness of neutrino masses is missing. 

If the SM and right-handed neutrinos reside in the bulk, the 5D Yukawa coupling 
X N is somewhat less constrained than the 4D coupling we were concerned with in 



the previous paragraph. The Dirac mass which enters eq. (|3.9| ) depends on the 
overlap of the neutrino zero modes and the Higgs 

m D = X N (H) = ^\% ) e-^H{y)^\y)fi N \y). (3.10) 

We assume that the Higgs profile has an exponential form which peaks at the TeV- 
brane, H(y) = H e 4k ^ nR \ This shape can be motivated by the equation of motion 
of a bulk scalar field [I|J]. Numerically it is equivalent to a delta function-like profile. 
Using the known mass of the W-boson, we can fix the amplitude H ~ Mpi in terms 
of the 5D weak gauge coupling. 

The zero mode of the SM neutrino has an exponential shape ( |2.6|) which for 
c v > 1/2 has only a small overlap with the Higgs at the TeV-brane. Sticking to 
the idea that all input parameters in the 5D theory should be of order unity in 
natural units, we assume that in the bulk lepton number is broken at M N ~ Mpi. 
Because of this large Majorana mass, the wave function of the lowest state of the 
right-handed neutrino is modified. Rather than acquiring a mass of order Mpi, this 
state gets localized at the TeV-brane where the warped geometry induces a typical 
KK mass for it, very similar to excited fermionic states. As a result, the overlap 
of the right-handed neutrino and the Higgs is large. These basic properties of the 
right-handed "zero mode" are fairly independent of the details of the KK reduction 
in the presence of a large bulk Majorana mass.0 They also do not depend on a 
possible 5D Dirac mass Cn for the right-handed neutrino. 

We now estimate the required value of X N to generate a 50 meV neutrino mass. 
As explained above, the Majorana mass of the zero mode Mjv is close to the KK scale 
Mkk- From comparison with electroweak precision data, Mkk cannot be smaller 
than about 10 TeV ||. By localizing /q close to the Planck-brane (c u ^> 1/2), 
rriD can be made very small. However, at the same time, the mass of the charged 
lepton, residing in an SU(2) multiplet with the neutrino, can become very small as 
well Since the tau is the heaviest lepton, it has to reside closest to the TeV- 

brane, and induces the largest rri£>. Assuming 5D lepton Yukawa couplings of order 
unity (i.e. ~ g^\ the 5D weak gauge coupling), the left-handed tau has to obey 
c < 0.57 (when the right-handed tau is delocalized, c =1/2) ||. If we take = 10 

k (6) 

,( 5 ) /„(5) On. m-5 



TeV, A^y = and c T = 0.57, we find mo = 7.9 GeV. A neutrino mass of 50 meV 



requires X N / ~ 9 x 10 5 . In accordance with the discussion in the previous 
paragraph, we have approximated the right-handed zero mode by the wave function 

4 If one allows for a suitably tuned profile of the Majorana mass in the bulk, the lowest state 
of the right-handed neutrino can be peaked around the Planck-brane. Its mass is still found to be 
bounded by the KK scale. 
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Figure 1: Wave functions of the v T and N zero modes, together with a sketch of the 
Higgs profile. 



of the first KK state of a fermion (|2.7f) without Majorana mass and = 1/2. In 
fig. [I] we present the associated wave functions of the v T and N zero modes, together 
with the Higgs profile. The bound on the Yukawa coupling only mildly depends on 
c^, increasing it from 1/2 to 1 would result in \ { ^/g {5) ~7x 1(T 5 . The error arising; 
from our approximation of the right-handed neutrino wave function is of the same 
level. 

(5) 

The bound on X N can be somewhat relaxed if we assume that the right-handed 
tau is localized at the TeV-brane (let's say with c = —1). Then we can move the 
left-handed tau and the tau neutrino closer to the Planck brane, c ~ 0.63, and find 
Aft/ 9^ ~ 6 x 10~ 4 . However, having the right-handed leptons localized towards 
the TeV-brane (c < 1/2) induces large deviations in the electroweak observables ||. 
The KK scale is required to be much higher then 10 TeV, and the gauge hierarchy 
problem gets reintroduced, albeit in a "mild" form. In the case of the muon the 
constraints on X N are an order of magnitude weaker. 

If one relies on supersymmetry to solve the gauge hierarchy problem, the Higgs 
fields can reside in the bulk. Then a different and more standard version of the see 
saw mechanism may be implemented: The right-handed neutrinos, having an inter- 
mediate mass, are confined to the Planck-brane. This translates into an intermediate 
see saw scale in the effective 4D theory. Since the Higgs is spread homogeneously in 
the bulk, there are no small overlaps anymore. Thus, an explanation of the fermion 
mass hierarchy is lost. An alternative way of implementing the see-saw mechanism 
is to assume a y dependent singlet Majorana mass. 

In conclusion, having the right-handed and sterile neutrinos in the bulk reduces 
the fine-tuning of the neutrino Yukawa coupling in the see-saw approach to some 
extent. However, a small number of order 10~ 4 is still needed to generate sub-eV 
neutrino masses. 
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4 Dimension Five Neutrino Masses 



In the 4D SM one does not have to worry about possible higher- dimensional oper- 
ators, since they are safely suppressed by powers of the huge Planck mass. Models 
with a low scale of gravity are very different in this respect. Naively one expects 
exotic non-renormalizable interactions suppressed by only a few TeV mass scale, 
unless they are forbidden by symmetries or multiplied by tiny coupling constants. 
In the warped SM the suppression scale for non-renormalizable operators can be 
anywhere between a few TeV and the Planck scale, depending on the localization of 
the fermions [||, 

As previously discussed in refs. @,|5[], Majorana masses for left-handed neutrinos 
are generated by the dimension-five operator 

J d'x J dy^g^HH iL C* jL = j d\ M^C*% (4.11) 

where are dimensionless couplings constants and C is the charge conjugation 
operator. The neutrino mass matrix reads 

< } = ^^e-^H\y)^\y)f^{y). (4.12) 

Again the most stringent constraint on / comes from the tau neutrino. Assuming 
the Yukawa coupling of the tau to be At 5 '' ~ and c(t r ) = 1/2, we found c(t l ) = 
0.567. From eq. ( [4.12|) we obtain a tau neutrino mass of m„ )T = I ■ 33 MeV 0. 



Bringing the neutrino mass down to 50 meV requires a very small value of / ~ 
10 -9 . Such a small coupling could originate from non-perturbative effects of gravity, 
especially if there is an extra dimension somewhat larger than Mpl [15]. However, 



it would be attractive to relax the bound on I within the framework of the effective 
theory we have considered so far. 

We can reduce the neutrino masses by shifting the left-handed leptons closer 
to the Planck-brane. To maintain the masses of the charged leptons, we either 
have to move the right-handed leptons closer to the TeV-brane or increase the 5D 
Yukawa couplings of the leptons A^ 5 -* . The first possibility is clearly disfavored 
by the electroweak precision data, as explained in the previous section. Keeping 
the positions of the right-handed leptons fixed at c = 1/2, the neutrino masses are 
roughly proportional to 1/(A^) 2 ||. If we restrict ourselves to = 10g^ 5 \ in 
order not to induce a large hierarchy in the 5D couplings, we find c(t~l) = 0.647 and 
m„ tT = I ■ 320 keV. 

The neutrino masses can be further decreased if we take the AdS curvature k 
to be smaller than the 5D Planck mass. Approximately, the relation m v oc (k/M 5 ) 2 
holds 0. The assumption k < M 5 is anyway necessary to derive the warped metric 
( |1 . 1|) from Einstein equations, neglecting terms with higher derivatives 0. In dis- 
cussions of the collider phenomenology of the warped model, the parameter range 
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M5/IOO < k < M5 has been considered ||16|| . Taking now the the very favorable 
case Ar 5 ' 1 = lOg^ 5 -* and k/M 5 = 0.01, we obtain m„ iT = 1-25 eV. A neutrino mass 
of 50 meV can therefore be generated by I = 0.002, which is only a moderately 
small number. However, one has to keep in mind that for kjM^ = 0.01 the radius 
kR = 8.63 is already of order 10 3 in units of the fundamental scale M 5 . 

Thus, the Majorana neutrino masses induced by the dimension five operator 
( |4.11| ) are quite large if we take the model parameters to be strictly of order unity. 
In ref. || we therefore imposed lepton number to eliminate ( [4. lip completely. Here 
we assume that one or a combination of the mechanisms discussed above is at work 
to bring down the neutrino masses to values which are consistent with experimental 
observations. In the following we investigate the consequences for neutrino mixing 
as well as for neutrinoless double beta decay. 



5 Neutrino Mixings 



In ref. H it was found that the dimension five neutrino masses ( |4.12| ) lead to small 
mixing angles for the neutrinos. This result arose from the different locations of the 
left-handed lepton flavors in the extra dimension in order to explain the hierarchy 
of charged lepton masses. In the following we demonstrate that large mixing an- 
gles naturally arise if the left-handed leptons have the same location in the extra 
dimension. This result is similar to the case of Dirac masses for neutrinos studied 
in ref. 
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We will also show that separating the electron neutrino somewhat from 
the muon and tau neutrinos helps to keep U e3 sufficiently small. 

To obtain the neutrino masses and mixings we diagonalize the mass matrix ([4.12 
with a unitary matrix = U U M^} C/J . The physical neutrino mixings 



77 M {u) U T 

u v Ivl di&g U v ■ 

U 



u}u v 



(5.13) 



also depend on the rotations of the left-handed charged leptons U\. The charged 
lepton masses arise from their couplings to the Higgs field in the same way as 
the Dirac neutrino masses of eq. ( |3.10|) , 

^ ^ ^e^H(y)f$(y)tf\y), 



M. 



(0 



2txR 



(5.14) 



where we have inserted the charged lepton wave functions and Yukawa couplings 
A; We determine Ui by diagonalizing M®(M^)l In general, the neutrino 

and charged lepton rotations are expected to show the same pattern on average 
since their positions in the extra dimension are tied together by the SU(2) gauge 
symmetry of the SM. Including CP violation, the mixing matrix U can be written 
as 

/ c 2 c 3 



U 



-C1S3 
S1S3 - 



- s 2 sic 3 e 
s 2 cic 3 e l 



id 



C1C3 - 
-S1C3 



C2S3 



s 2 c 1 s 3 e" 



s 2 e 

C 2 Sl 

c 2 ci 



\ 



(5.15) 
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where Sj and q are the sines and cosines of d{. Like in the CKM matrix, there 
is a single complex phase S which induces CP violation in the lepton sector. The 
neutrino mass matrix contains two additional Majorana phases which, however, do 
not show up in the mixing matrix U. 

The atmospheric neutrino data imply |jT7 



1 • l(T 3 eV 2 < Aml tm < 5 • l(T 3 eV 2 , sin 2 29 1 > 0.85. (5.16) 
There are several solutions to the solar neutrino anomaly [T^ 



Am s 2 ol [eV 2 ] sin 2 26 3 

LMA 2 • KT 5 - 4 • 10~ 4 0.3-0.93 

SMA 4 • 10~ 6 - 9 ■ 10- 6 0.0008- 0.008 (5.17) 

LOW 6 • 10~ 8 - 2 • 10~ 7 0.89 - 1 

VAC ~ 10- 10 0.7-0.95 

The SMA solution gives only a poor fit to the data. The CHOOZ reactor experiment 
constrains \U e 3\ 2 = s\ to be at most a few percent [19]. Nothing is known about the 
CP violating phase 5. 

Our solution to the neutrino puzzle follows the "neutrino mass anarchy" mod- 
els ||20|| . Mass matrices with randomly chosen entries of order unity have a large 
probability to fit the neutrino data. Thus, we parametrize the coefficients in (|4.12 ) 
by 

Zy = A.f fi . (5.18) 

In contrast to the 4D realization of "neutrino mass anarchy", in our scenario the 
smallness of the neutrino masses is not attributed to a very small overall magnitude 
A of the couplings. Rather, the neutrinos are localized close to the Planck-brane, 
where the dimension five operator Qi.llj) is suppressed. The supposedly fundamental 
theory responsible for the effective interaction ( 4.11 ) is represented through the order 



unity coefficients Uj. To incorporate the charged lepton rotations we also generate 
charged lepton mass matrices with random Yukawa couplings of order unity. We 
choose suitable lepton locations to reproduce the measured lepton masses. 

We consider the very favorable case of \ { r ] = 10g (5) and k/M 5 = 0.01. We take 
the position of the right-handed tau to be c t> r = 1/2, which means its wave function 
is delocalized in the extra dimension. To reproduce the tau mass we find the position 
of the left-handed tau to be c Ti l = 0.72. Since different positions for the left-handed 
leptons tend to produce small neutrino mixing angles @, 0, we take 

c e , L = c^ L = c TtL = 0.72. (5.19) 

To accommodate the muon and electron masses we use c^r = 0.64 and c e ^ — 0.85. 
Taking the order unity coefficients to be homogeneously distributed 1/2 < < 2 
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Figure 2: Distribution of \U e3 \ 2 for the case c e ^ = c^l = c Tj l and 1/2 < \l^\ < 2. 



with random phases from to 2tt, we find the most favorable value for the overall 
scale to be A = 1.9- 10~ 3 . We randomly generate parameter sets for Z^-, calculate the 
neutrino mass matrix from (|4.12|) and compute the neutrino masses and mixings. 
For the charged lepton rotations we take random Yukawa couplings 10p^/\/2 < 
I [^i ]ij I < 10(?^ 5 ^ ■ y/2 with phases from to 2n. The charged lepton masses and left- 
handed rotations are then calculated from the mass matrix eq. (|5.14 ). We require 
the computed lepton masses to agree up to a factor of 3/2 with the measured values, 
which holds for about 50% of the random sets of Aj . These sets are the starting 
point for the investigation of the neutrino properties. 

Focusing on the LMA solution of the solar neutrino anomaly, which turns out to 
be clearly favored, we find the following picture. From our neutrino parameter sets 
about 70 % reproduce Am 2 tm ( |5.16|) . Imposing in addition the constraint from Am^, 
( |5.17| ), we are left with about 28 % of the parameter sets. The solar and atmospheric 
mixings angles bring this number down to about 6 %, which is still a considerable 
fraction given the number of constraints. The most stringent constraint turns out 
to come from the CHOOZ experiment. Adding the requirement \U e s\ 2 < 0.05, the 
fraction of viable parameter sets shrinks to about 0.7 %. This result is clearly 
related to (\U e 3\ 2 ) — 0.22, which is considerably above the experimentally favored 
value. Here we use the logarithmic average of a quantity X 



(X) = exp 



N 



N 



(5.20) 



If not stated otherwise, we average over the parameter sets which reproduce the 
correct Am 2 's. In fig. ^|we display the distribution of \U e 3\ 2 found in our statistical 
analysis. We included only parameter sets that reproduce the correct Am 2 's. Small 
values of \U e 3\ 2 are somewhat favored. The corresponding distributions for the solar 
and atmospheric mixing angles are peaked at maximal mixings. Averaging only 
over the sets which pass all experimental constraints, we find (\U e 3\ 2 ) = 0.022. This 
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^ m atm,sol 


+ Sill 2 26 , atm ,sol 


+ \U e3 \ 2 < 0.05 


LMA 


44.4 (28.1) 


5.8 (5.9) 


5.0 (0.7) 


SMA 


1.3 (0.04) 


0.3 (< 0.001) 


0.3 (< 0.001) 


LOW 


0.008 (< 0.001) 


0.002 (< 0.001) 


0.002 (< 0.001) 



Table 1: Probability in percent that a randomly generated set of coefficients 1/2 < 
\lij\ < 2 satisfies the constraints from Am 2 tmsol (first column) and sin 2 26 , atmiSO i 
(second column) and \U e ^\ 2 < 0.05 (third column). The results are given for the 

Case C etL > C^ L = C T:L (c e:L = C^ L = C TjL ) . 



means that \U e 3\ 2 is probably close to the experimental bound and can very likely be 
tested at future neutrino experiments, such as MINOS, which is sensitive down to 
| f7e3 1 2 ~ 0.0025 [21]. Since we typically start with large phases in the neutrino and 
charged lepton mass matrices, the complex phase 5 in the mixing matrix U is found 
to be most likely of order unity. If we ignore the charged lepton rotations, i.e. take 
U = U v , the fraction of parameter sets satisfying all constraints rises somewhat from 
0.7 % to 1.7 %. The reason is that cancellations amongst the large mixing angles in 
U\ and U u make it slightly more difficult to reproduce large solar and atmospheric 
mixings. 

These results depend only mildly on the interval chosen for the coefficients Uj. 
If we narrow the range, let's say to 0.8 < < 1.7, the fit for the Am 2 's, the solar 
and atmospheric mixings angles improves remains unchanged, while the CHOOZ 
constraint becomes slightly harder to satisfy. We find that now about 0.6 % of the 
parameter sets meet the constraints from eqs. ([5.171) and (|5.16|) and \U e3 \ 2 < 0.05. 



The SMA, LOW and VAC solutions to the solar neutrino anomaly can be realized 
only with severe fine-tuning. The required small values of Am 2 ol are very unlikely 
to be produced through accidental cancellations in the neutrino mass matrix. For 
the SMA we find that for the most favorable value of A = 1.6 • 10~ 3 and 1/2 < 
< 2, only about 0.04 % of the parameter sets reproduce the correct Am 2 's. 
Including the constraints from the solar and atmospheric mixing angles and the 
CHOOZ experiment brings the fraction down to less than 10~ 5 . The LOW and 
VAC solutions are even more fine-tuned. Under the given assumptions the LMA 
case is therefore selected as the only realistic solution to the solar neutrino problem. 
These results, which are summarized in table [II agree nicely with the 4D model 
studied in ref. Note that for the LMA solution the results remain stable if 



we use real instead of complex mass matrixes. However, the SMA, LOW and VAC 
solutions are more difficult to accommodate with complex mass matrices since the 
required cancellations to obtain a small Am 2 ol become even more unlikely. For 
instance, with real mass matrices, the fraction of parameter sets which reproduce 
the SMA mass squared differences increases from 0.04 % to 0.7 %. 

So far the constraint from the CHOOZ experiment is the most stringent one. 
Its inclusion reduces the probability that a parameter set realizes the LMA solution 
from about 6 % to less than one percent. With all entries in the neutrino mass matrix 
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Figure 3: Distribution of \U e ^\ 2 for the case c e ^ > c^i = c Tj l and 1/2 < \l^\ < 2. 

being of similar magnitude, the ensuing mixing angles are typically large. The fit 
to the neutrino data improves considerably if the electron neutrino is somewhat 
separated from the other two neutrino species. Shifting the electron neutrino closer 
to the Planck-brane induces small elements in the neutrino mass matrix. As a result, 
small values of Am^, and \U e ^\ 2 become more probable. The neutrino mass matrix 
Q4.12Q acquires the following structure 





( ' 2 


e e\ 




e 


1 1 




V e 


1 W 



where e « f^J (ttR)/ foli^R) ~ exp(— (c Si l — c^^nkR). The charged lepton mass 
matrix does not follow the pattern (|5.21| ) because of the different locations of left- 
and right-handed leptons. Still U\ and U u 

We keep the muon and tau neutrinos at the previous locations c^l = c r ,L — 
0.72. A separation of these two would make the fit to the atmospheric neutrino 
data more difficult since the corresponding mixing is reduced. (A small separation, 
e.g. \c^l — Ct,l\ < 0.02 could, however, be tolerated.) For the electron neutrino the 
most favorable choice turns out to be c e ^ = 0.79 leading to e = 0.15. This value 
is of the order of y / m jU /m r , and also of the Cabbibo angle. However, in our model 
there is no relation to either of these quantities. The positions of the right-handed 
leptons are found to be c 6i # = 0.79, c^r = 0.62 and c Tj ^ = 0.49. 

Let us again focus on the LMA solution. Assuming 1/2 < < 2 for the 
order unity coefficients, we find the best value A = 2.6 ■ 10~ 3 . The correct values of 
the Am 2 's are fitted by 44 % of the parameter sets. Taking into account also the 
constraints from the solar and atmospheric mixing angles, this fraction shrinks to 
5.8 %. This reduction is mostly due to the solar mixing angle which is suppressed 
if e is small. Most important, however, the CHOOZ constraint \U e 3\ 2 < 0.05 is now 
satisfied almost automatically, and we are finally left with 5.0 % of the parameter 
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sets. Compared to the case of e = 1, the probability for a parameter set to satisfy all 
observational constraints is enhanced by almost an order of magnitude. It turns out 
that including phases in the mass matrices is very important to obtain this result. 
With real valued mass matrices there tend to be large cancellations between the U u 
and U\ contributions to the solar and atmospheric mixing angles. Large mixings 
are then very difficult to realize. If we consider the interval for the coefficients to 
be 0.8 < < 1.7, the fit becomes slightly better because the solar mixing angle 
gets somewhat enhanced. We find that 5.5 % of the parameter sets pass all the 
constraints. 

In fig. |3| we display the distribution of \U e s\ 2 for the case 1/2 < < 2, where 
we again include only the parameter sets which reproduce the correct Am 2 's. As 
expected, the distribution is now peaked at small values of |?7 e 3| 2 - Averaging only 
over parameter sets which satisfy all constraints, we find a moderately small value 
(\U e 2,\ 2 ) — 0.010. As a consequence, the next generation neutrino experiments still 
has a good chance to detect a non- vanishing \U e z\ 2 |]21] . 



Taking even smaller values of e, it is possible to implement the SMA solution to 
the solar neutrino problem. The most favorable choice of parameters we find to be 
c 6j l = 0.96 corresponding to e = 0.0020 and A = 2.6 • 10~ 3 . The correct values of 
Am 2 are reproduced by 1 % of the sets of coefficients 1/2 < < 2. The constraints 
from the solar and atmospheric mixing angles reduce this amount to about 0.3 %. 
The CHOOZ constraint is always satisfied, and we obtain (|t/ e 3| 2 ) = 2.7- 10~ 6 which 
is much too small to be measurable. The LOW solution is very difficult to implement 
because of the small solar Am 2 . We take c e ^ = 0.86 corresponding to e = 0.0026, 
and A = 2.4 • 10~ 3 . The constraints from the Am 2 's are satisfied by only about 0.01 
% of the parameter sets. All the neutrino data are reproduced by a fraction smaller 
then 10 -5 . The CHOOZ constraint is again satisfied automatically. Small Am 2 's 
could be accommodated by moving the electron (and muon) neutrino closer to the 
Planck-brane. Then it becomes however more difficult to obtain large mixings. The 
VAC solution is even more difficult to realize because of the very small value of 
Am 2 ol . Thus, the LMA solution is by far the most favored scenario. A collection of 
our results is given table [1| 

Neutrino mass matrices of the type ( |5.21| ) have previously been considered in 



refs. [22, 23]. There the small quantity e was attributed to a Froggatt-Nielsen mech- 



anism 



24] . Our findings agree very well with the results of these studies. 



6 Corrections from KK Neutrinos 

So far we have only considered the neutrino zero modes. But eq. ( |4.11| ) also induces 
mixings between the zero modes and the vector-like excited states. In the following 
we show that this effect does not modify the conclusions we have reached above. 
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The general neutrino mass matrix takes the symmetric form 



/ 



Mr 



(o) 



V 



R ' 



m 



m 



(0,0) 
L 

(1,0) 



m 



in 



(0,1) 
L 

(1,1) 



rriKK l Tti 





mKK,i 
(i,i) 



\ 



V 



V 



V 



(1) 

L 

(1) 
R 



(6.22) 



" • / V : / 



Here m 



(0,0) 



are the Majorana masses for the zero modes from eq. ( [4.121) . 



The Majorana masses ml and m^ 31 are obtained from eq. ( |4.12| ) by inserting 



the wave functions of the relevant states. m^Kj are the KK masses of the excited 
states. Since the KK states have a greater overlap with the Higgs, the corresponding 
Majorana masses are larger than those involving zero modes. We consider again the 
case k = 0.01M 5 and At = 10g^ implying c t> l = 0.69. To estimate the magnitude 
of the effect we restrict ourselves to a single neutrino flavor. For the first excited 



state we find m 



(o,i) 



I ■ 8.8 keV, m 



(i,i) 



I ■ 3.5 MeV and m 



(i,i) 

R 



I ■ 20 keV. For 



higher excited states these masses approximately do not change. If we would have 
used a delta function-like Higgs profile, would be zero since the wave functions 
of are odd and vanish at the boundaries. Truncating the KK tower at the first 
excitation, we find for the lowest state (i.e. the zero mode) 



,(°) 



+ 



m 



(i,i)™(°>i) 

R 



m 



^ - 



m 



(0,1) 



m 



KKA 



rriKKA 



4" 



(6.23) 



v 



(0) 



+ 1.5- 10 



- 18 / 2 z/i 1} -8-10- 



10 ^ (1) 



R 



In the second line we used numbers from the previous example. These tiny admix- 
tures have only a negligible input on the properties of the zero modes. The tiny 
correction to the mass of the zero mode, for instance, is given by 



5m,, 



in 



(1,1)^(0,1)^2 
R 



[m 



(6.24) 



1.3-10 



KKA. 
14 Z 3 



eV. 



The contributions of the higher KK states are suppressed by larger KK masses and 
will affect our result by at most a factor of order unity. In a similar way to ( |6.22| ) the 
charged lepton masses (|5.14j) also couple zero modes and KK states. Here too the 
large KK masses suppress corrections to charged lepton masses and mixings from 
the excited states. Thus, the KK states can be safely neglected in the discussion of 
neutrino mixings. 



7 Dirac vs. Majorana Neutrinos 

We finally discuss some experimental signatures to distinguish the presented scenario 
of Majorana neutrinos from the Dirac neutrino scheme discussed in refs. |12,|i~3| . 
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Figure 4: Fraction of parameter sets p which reproduce the Am 2 's as a function of 
the electron neutrino Majorana mass m ee for cl = 0.72 (in logarithmic scale). 



The most direct evidence for a Majorana mass of neutrinos would be the discov- 
ery of neutrinoless double beta decay (0z//3/3). The non-observation of this process in 
the Heidelberg-Moscow experiment implies an upper bound on \ = m ee < 0.35 



eV ||25|| . There are plans to bring this limit down to about 0.01 eV ||26|| . If all neu- 
trinos are at the same position in the extra dimension, we expect m ee ~ 0.02 eV, 
which is far below the current experimental sensitivity but within reach of future 
experiments. This result does not depend on what solution to the solar neutrino 
anomaly is realized. Fig. [| illustrates the unnaturalness of a large Majorana mass 
in our scenario. As a function of m ee we show in a logarithmic scale the fraction of 
parameter sets p which reproduce the Am 2 's required for the solar LMA solution 
and for the atmospheric neutrino anomaly. (Including also the constraints from the 
mixing angles would reduce the number of viable parameter sets by another order of 
magnitude.) We assume all neutrinos to be at the same location with cl = 0.72, and 
take 1/2 < < 2. To tune m ee we set In = 1 and vary A. For small values of m ee 
(i.e. small A) it becomes impossible to accommodate the atmospheric Am 2 , leading 
to p = 0. For large Majorana masses p becomes exponentially small since it is highly 
improbable to accidentally generate Am 2 -C m 2 . The peak around m ee ~ 0.015 eV 
corresponds to the fraction of 28 % quoted in table [|. For larger ranges of lij the 
distribution is spread out further. 

Once the electron neutrino is localized closer towards the Planck-brane Qv(3(3 
becomes drastically suppressed proportional to e 2 . For the clearly favored LMA sce- 
nario we find m ee ~ 0.001 eV, which is too small to be detected in the near future. 
For the SMA, LOW and VAC solutions m ee is even smaller. It is therefore ques- 
tionable if our model induces 0z//3/3 at a detectable level. In any case the Majorana 



masses we find are far below the recently claimed evidence of about 0.4 eV |27|. 

Tritium endpoint (beta decay) experiments are sensitive to neutrino masses of 
Majorana and Dirac type. The MAINZ collaboration plans to look for neutrino 
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masses down to about 0.3 eV p8 |, which is , however, still an order of magnitude 
larger than the typical neutrino masses in the scenario considered. 

In refs. |5|, [13|] we imposed lepton number symmetry to make the proton stable 
and eliminate Majorana neutrino masses from non-renormalizable operators. If in 
the current framework we want to keep Majorana masses but still forbid proton 
decay, lepton parity is an attractive possibility. It still would allow baryon number 
violating processes, such as neutron anti-neutron oscillations, which might be close 



to the observational bound [jnj. Without some symmetry dimension six operators 
inducing proton decay have to be multiplied by small couplings of order 10~ 8 . The 
proton lifetime then should not be too far above the present experimental constraint. 

The scenario with Dirac fermions could manifest itself through lepton flavor 
violating decays. For instance, the rate for fi — > cy is considerably enhanced by the 
presence of sterile neutrino KK states, which spoil the GIM cancellation of the SM. 
If the SM neutrinos are confined to the TeV-brane, the branching ratio for /i — > 
is above the experimental limit of about 10 -11 , unless the KK scale is above 25 TeV 
2Pj . In the case of bulk SM neutrinos fi — > ey does not constrain the KK scale 



anymore, but the branching ratio may still be close to the experimental bound. 



8 Conclusions 

Neutrino masses generated by dimension five interactions in the SM are of magnitude 
< 10~ 5 eV, and consequently the atmospheric neutrino anomaly remains problem- 
atic. In contrast, such masses are significantly larger in the 5D warped SM, and we 
have shown that an explanation of the solar and atmospheric neutrino anomalies can 
be realized, based on the idea of neutrino mass anarchy. We also provide estimates 
for the mixing angle U e 3 as well as the parameter m ee that appears in the amplitude 
for neutrinoless double beta decay. No new particles such as SM singlet neutrinos 
are needed to implement the scenario. Indeed, because of the warped geometry, the 
usual see saw mechanism will not work in its simplest form, without either giving 
up the resolution of the gauge hierarchy problem, or invoking supersymmetry. 
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